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Abstract: We show how the scheme- and scale-dependence of the short- distance 
operator product expansion with four-quark operators can be correctly accounted for 
in the framework of the 1 /A^c-expansion once the hadronization of two-quark currents 
has been fixed. We show formulas explicitly in the case of the S/^-parameter. We 
then use them with our earlier estimates of the long-distance effects. We compare 
Chiral Perturbation Theory at Leading- and Next-to-Leading-Order with the ENJL 
model results in the chiral limit. Good matching between the long- and short-distance 
regimes is obtained and our final value for the physical scheme independent Bk is 
0.77 ±0.05 ±0.05. 
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1. Introduction 



Weak non-leptonic decays and mixings are one of our few windows on the CP- 
violating sector of the Standard Model and their calculation is thus important. A 
set of reviews summarizing the present status in the kaon area can be found in the 
various talks at the workshops at Orsay and Chicago[l, 2]. 

The large difference in mass between the PF-boson and the kaon presents an 
additional difficulty since logarithms of this ratio need to be summed to all orders. 
This can be done using the Operator Product Expansion (OPE) and is by now 
standard. The remaining problem is the calculation of the matrix elements of these 
operators at some low scale. 

That the l/A^c-^xpansion would be useful in this regards was first suggested by 
Bardeen, Buras, and Gerard[3, 4] and has been reviewed recently by Bardeen[5]. 
There one can find most of the references to previous work and applications of this 
non-perturbative technique. 

Since the original work [3] there are have been some improvements. This has 
centered on identifying more correctly the scales in short- and long-distance[6, 7] 
and their matching. A more precise method to perform this identification was pro- 
posed in [8], the X-boson method, and shown there to provide a correct matching 
with the renormalization group evolution at one-loop. At ncxt-to-lcading (NLO) in 
the renormalization group running several other problems appear. The operators 
become scheme dependent and also dependent on various other choices like the one 
of evanescent operators. This is discussed in the review by Buras [9] and also in the 
original OPE at NLO papers [10, 11]. 

In this paper we show how scheme and scale dependence can consistently be 
treated within the 1/A^c-expansion technique as argued in [4, 5, 12, 13]. We will 
use the method suggested in [8] as already used in [12, 13]. In particular we will 
show explicit expressions for the current x current AS — 2 four-quark operator 
in the Standard Model. In practice this means tracking the scheme dependence 
consistently during the entire calculation. 

We are concerned here with the short-distance matching between the definition 
of the weak operators done in perturbative QCD using a particular regularization 
scheme and the actual calculation of weak matrix elements at a given order in the 
1/Nc expansion. This can be done purely within perturbation theory and is unam- 
biguous once the short-distance schemes used in perturbative QCD are specified. We 
will also show explicitly that it docs not depend on the infrared regulators. We then 
use this result together with the matrix element obtained earlier in [12] to obtain 
complete scheme-independent results for the parameter within the 1/A^c expan- 
sion. We also discuss the long-distance short-distance matching for this quantity. 
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2. The Flavour Changing Effective Action and the X-Boson 
Method 

In light hadrons flavour changing decays there appear two very different scales; 
namely, the hadronic scale and the weak scale, which makes necessary the use of 
effective held theory analyses. This is done with the help of the OPE which allows 
to separate short-distance from long-distance physics. The OPE analysis is done at 
some scale below the W mass and this separation introduces some short- distance scale 
and regulator scheme dependence which has to cancel in the final physical amplitude. 
We will sec how this occurs explicitly in the l/N^. technique. 

The process of obtaining the relevant effective action continues by integrating 
out the heavy degrees of freedom up to some scale below the charm quark mass 
and setting the appropriate matching conditions at heavy particle thresholds [9]. 
This can be done using QCD perturbation theory. In this process, large logarithms 
{as{i') ln{Mw/i')) need to be resummed to all orders in perturbation theory. This 
is done with the help of Renormalization Group (RG) techniques which at Ncxt-to- 
Leading-Log-Order (NLLO) allow to rcsum up to as(z/) {as{i')hi{M\Y/i/))'' . In fact, 
only from NLLO can a full scheme dependence study be done. 

The whole process described schematically above is by now standard and has 
been brought up to NLLO in [10, 11] from the earlier Leading-Log-Order (LLO) 
results[14]. For comprehensive reviews where complete details can be found, see [9]. 

Thus, at some scale v below the charm quark mass, we are left with the effective 
field theory action 

rA5=a = -CA5=aE<^^W [ V Qiiy)+^-<^- (2-1) 

to describe hadronic fiavour changing processes in the Standard Model. Here Qi{x) 
is the set of four-quark operators made out of the three light-quark fields. They 
change fiavour in a units and are defined perturbatively within QCD using some 
regularization. In particular, the NLLO calculations mentioned above have been 
done both in the 't Hooft-Veltman (HV) scheme (MS subtraction and non-anti- 
commuting 75 in D 7^ 4) and in the Naive Dimensional Regularization (NDR) scheme 
(MS subtraction and anti-commuting 75 in 7^ 4). 

In order to reach a physical process we now need to calculate matrix elements 
of the effective action (2.1) between the physical states. In this calculation, all 
dependence on the different schemes used and on the scale u should disappear. This 
necessarily involves long-distances which cannot be treated analytically within QCD. 
One option is to simply give up here and refer to lattice QCD calculations. However, 
these have at present many problems with this type of matrix elements. They are 
complicated quantities to treat on the lattice and chiral symmetry effects are very 
important. For a recent review on lattice QCD weak hadronic matrix elements see 
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e.g. [15]. Here we would like to show how to connect the effective action in a well 
defined way with other approaches to low-energy QCD. 

To identify a four-quark operator in any model or approximation to QCD like e.g. 
Chiral Perturbation Theory (CHPT), is very difficult. On the other hand, all these 
approximations and models are normally tuned to reproduce various experimentally 
observed properties involving currents and/or made to specifically match perturba- 
tive QCD predictions for Green functions involving currents. As a rule therefore, 
identification of two-quark currents is possible. So a first step is to replace rAS=a by 
an equivalent effective action that only involves currents. 

We thus introduce an action of X-bosons with masses Mx coupling to quark 
currents, whose OPE reproduces Eq. (2.1). This requirement fixes the couplings 
of the X-bosons. This matching can be done at a perturbative level in QCD with 
external states consisting of quarks and gluons as long as the scale u is high enough. 
At this matching the scale and scheme-dependence of Eq. (2.1) has disappeared and, 
as we show below, the X-boson couplings are independent of the precise infrared 
regulator chosen here. 

We now split the integral over X-boson momenta into two parts. The high 
energy part needs to have the momentum flow back through perturbative quarks 
and gluons leaving an operator behind that only needs to be evaluated to leading 
order in l/Nc- The low energy part needs to be evaluated using the 1/Nc method 
to NLO in l/Nc- Putting the two parts together we can see that all dependence 
on the X-boson masses has dropped out and the final result can be written in a 
way that only involves the coefficients of the operators in Eq. (2.1) with the scheme 
dependence removed correctly. In the next section, we will show this procedure 
explicitly for the AS" = 2 four-quark operator. The same procedure can be worked 
out for the AS — 1 effective action but is much more cumbersome[16]. 

3. The AS' = 2 Case 

We present a complete analysis for the A^* = 2 case. The generalization to other 
fiavour changing processes is straightforward. The Standard Model effective action 
for AS — 2 transitions is 

rA5=2 = -CaS=2C{u) Jd^y QAS=2{y)+h.C. (3.1) 

with 

Qas=2{x) = AL^{x)L^{x) ■ 2L^ = [s^ (1 - Ts) d] {x) . (3.2) 

The operator Qas=2 transforms under SU(3)l x SU(3)r rotations as a 27l x 1r. 
The normalization factor 

Cas=2 = ^Hml ml, M^, Vckm) (3.3) 
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is a known function of the integrated out heavy particles masses and Cabibbo- 
Kobayashi-Maskawa matrix elements and the Wilson coefficient C(i/) is known to 
NLL0[11]. 

The AS" = 2 K'^-K^ matrix element 

{K^\K') = -iCas=2C{u){K^\ J d^y QAS=2{y)\K') (3.4) 

governs the short-distance contribution to the Kl-Ks mixing. We describe now how 
to calculate it in the 1/Nc expansion [3, 8, 12, 17]. 

3.1 Transition to X-boson Effective Theory 

First of all, wc need to set the matching conditions between the weak operator Qas=2 
defined by (3.2) and the short-distance part of the effective field theory of X-bosons 
used in the 1/Nc technique. The effective action of the latter is 

Tx = 2^A5=2(//c, Mx, • • •) / X^{y) L^{y) + h.c. (3.5) 

where gAS=2{t^c, Mx, ■ ■ ■) is an effective coupling. In addition we add kinetic and 
mass terms for the X-bosons. Evaluating matrix elements of rAs=2 needs to be done 
in the scheme in which it is defined. We can freely choose the scheme in which we 
should treat Fx- Here we choose an Euclidean cut-off regulator with subtraction to 
define the strong coupling ag at a scale /ic- This has certain advantages for the next 
subsection. 

We now determine the coupling gAS=2 by requiring that the matrix elements of 
rA5=2 are the same as those of Fx to leading order in 1/Mx- We require at some 
perturbative scale 

{s{q,)d{q^)\e'^-^=-\s(qMqs)) = {s{q,)d{q^)\e^^-\s(q2)d{qs)) + C»(l/M^) ■ (3.6) 

Where we have taken the matrix elements between an incoming strange anti-quark 
with momentum q2 and an incoming down quark with momentum ^3 and an outgoing 
strange quark with momentum qi and an outgoing down anti-quark with momentum 
q4. We regulate infrared divergences by having — 5? > 0. In order to have expressions 
of manageable length we require —qf <^ jg^ • qj^i\ <^ //^,i/^,M|- but we have kept 
all momenta different in order to show the full cancellation and remove any possible 
ambiguities in applying the equations of motion. 

Both sides in (3.6) are by themselves scale, scheme, and regulator independent 
at at given order. This is to order a^(z^) if we use NLLO results. They do depend 
on the infrared regulators used, like the quark momenta, and the gluon gauge but 
these are the same in both sides and cancel in the matching. The gauge dependence 
we discuss in Section 3.3. 

The left-hand side of (3.6) evaluated at NLLO using the diagrams of Figure 1 

is 
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(b) (c) (d) 

Figure 1: The type of diagrams contributing to the matrix element of Ti\s=2- The 
circle denotes the current-current operator. Colours are connected on the same side of the 
horizontal line. The curly line is a gluon and full lines are quarks, (a) Wave Function 
Renormalization. (b) Vertex Diagrams. (c)(d) Box Diagrams. 



with 



1 + q,, qs, q,)^] S^+(l + F{q,, q^, -q,, -q^)^] S2 



7i /^2gi.g2 



ri 



(3.7) 



where i/ is a scale where perturbative QCD can be used. The function 



1 N"^ -1 

F{qi, q2, gs, 54) = B{qi, gs, qs, Qi) - J^^^l^^ ^2, -^4, -Qz) + V{qi, q2, qs, Qi) ] 

(3.9) 

collects finite terms from the box (B) and vertex diagrams (V). 



S, = [si7m(1 - 75)^3] [^27^^(1 - l5)dA] 
S2 = [si7m(1 - 75)c?4] [527^^(1 - 75)4] 



(3.10) 



are tree level matrix elements. The colour indices are summed inside brackets. The 
\E'j field destroys a quark with momentum g^. In the NDR scheme, using the 
Feynman gauge for the gluon, the one-loop anomalous dimensions is 



and 



7i 



^NDR 



2 V Nr 



4 V Nr 



The one-loop anomalous dimensions 71 is scheme-independent and 

rr = --(i-- 



(3.11) 
(3.12) 

(3.13) 



if we use the 't Hooft-Veltman scheme instead. In both cases we have used the same 
scheme for the evanescent operators as the one in [11, 18]. 

The terms with F{qi, q2, qs, q^) collect the infrared dependence on quark masses, 
external quark momenta, ■ ■ ■ as well as scheme-independent constants. If we had 
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used a gluon mass as infrared regulator that would have been present there as well. 
The separation of the constant terms between ri and F{qi, q2, qs, q^) is arbitrary. 
The explicit expression for the box function B{qi, q2,q3, q^) is 



B{qi,q2,q3,qA) 




+ In 



-2qi ■ q2 



In 



-2gi ■ q2 



- In 



gi . 



hi 



-2^2 • q\ 



ql 



+ 0{qt/{q^ ■ q,^j)) 



(3.14) 



and for the vertex function V{qi, q2, ^3, 54) is 



V{qi,q2,q3,q4) = 



2+l|! + 21n(z^]lnp^-^^ 



+ 2in(^%:^)iM^^)-^- 



?2 



?4 




• '13 



_3 / -2^2 -94 ^ 3 / -2^2 • g4 
2 I gi j 2 I g| ^ 



+ 0{q^/{qj ■ qk^,)) . 



(3.15) 



In the HV scheme one gets the same expressions for these two functions. 
The Wilson coefficient C{v) at NNLO can be written as 



C(^)= 1 + 



TT 



72 I32ll 



71/^1 



(3.16) 



with [11] 



7?°" = 

















57 / AA2^ 



-17 + 4(n/-3) + 



57 



NT 
9 



16A 



(3.17) 



with Pi — —9/2 and /32 = —8. The Feynman gauge for the gluon has been used to 
obtain the results above [11]. The scheme of evanescent operators we used is the one 
of [11, 18]. 

Equation (3.7) is now scheme-independent and scale-independent to order a^, 
the difference in running via Eq. (3.17) is compensated by the difference between 
(3.12) and (3.13). This was shown with simpler states also in [18]. We can now 
check this by plotting —Cd/Cas=2 versus u for the two schemes. The difference 
and the variation with u is an indication of the neglected corrections. This is 
shown in Figure 2 where we plot the one- loop coefficient C{u) and the two- loop one 
for both the NDR and the HV schemes. They differ considerably and are rather 
i/-dependent. Now —Cd/Cas=2 of Eq. (3.8) is more stable with u and has a smaller 
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scheme-dependence left. Both are an indication of the size of the corrections. This 
can be improved systematically by doing the matching to order and the running 
to three- loops and so on. We have chosen 2qi ■ q2 = I GeV^ and ^^(M^) = 0.334 [22] 
as input. 
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Figure 2: The short-distance coefficients. 

For the left-hand side of (3.6) we need to specify now which is the effective field 
theory we use in the l/Nc calculation. We introduced a fictitious [4, 8] boson 
which reproduces the physics of the Qas=2{x) weak operator below fic ~ i^, with 
coupling given in Tx of Eq. (3.5). 

For the right-hand side of (3.6) at NLLO we get 



iC, 



c 



I + F{qi,q2,q3,q4 



c) 



Si+ [1 + F{qi, q2, -q4, -qs. 



as(/i 



c) 



TC 



+0 



S2 

(3.18) 



using the diagrams of Figure 3, with 



C, 



c 



'9as=2 
Ml 



1 + 



«s(/i 



C 



TT 



.2 V ^1 . 



+ ri 



(3.19) 
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Figure 3: The type of diagrams contributing to the matrix element of X-boson exchange. 
The wiggly line is the X-boson, the curly line is a gluon and full lines are quarks, (a) Wave 
Function Renormalization. (b) Vertex Diagrams. (c)(d) Box Diagrams. 



We have used the conditions in the X-boson effective theory Mx ^ qf and fiQ ^ 
—qf- For the scheme dependent constant we get 



The box diagrams are finite in this case, the only divergent ones are the vertex and 
wave-function renormalization diagrams, (a) and (b) in Figure 3. 

The expression in (3.18) is also scale- and scheme-independent at order a^. The 
constant fi is scheme^ dependent but it cancels against the same dependence in the 
constant gAs=2- 

About the QCD regularization used, notice that there are no logarithms depend- 
ing on the scale fic- This is not trivial, in principle these can also appear. This is a 
consequence of the fact that the anomalous dimension vanishes for two-quark vector 
and axial-vector currents and that box diagrams are finite. Two-quark currents can 
be identified always in the low-energy approximation to QCD unambiguously as said 
before. And for box diagrams we can use D = 4 and therefore their contribution is 
75 and evanescent operator schemes independent. 

The infrared regulator scheme dependence is , as it should be, precisely the same 
in Eq. (3.7) and (3.18), precisely the same function -F(gi, 52, ?3, ^4) appears. 

The X boson mass acts here like an ultraviolet regulator. The dependence has 
to disappear in the final physical amplitude. We will see how this happens in the 
matching between long- and short- distances. 

So from the matching condition in (3.6) we obtain for the X-boson couplings: 



^Hcre scheme means the precise definition of the X-boson coupUngs and QCD regularization 
used. 




(3.20) 



9ls=2if^c, Mx, ■■■)= MICas=2C{v) 1 + 
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and as{i')—asifJ'c) is also of order a^. Wc also plotted —g\s=2/ {^x^ ^3=2) in Figure 
2 to show the effect of fi for the HV and NDR schemes. Again we see that the result 
is stable with v. We used Mx = 1 GeV for definiteness. 

With this first matching we have obtained an expression which is scale and 
scheme independent to NLLO order. In addition in the X-boson theory the only 
part that needs regularization is the X-boson-quark vertex itself and we have not 
used any 1/Nc argument until now. 

3.2 Long-Distance Short-Distance Matching and 1/Nc 

With the value of gAS=2 set in (3.21) we can calculate the weak matrix element in 
the 1/Nc expansion within the X-boson exchange effective theory. 

We follow the same technique we used for calculating the electromagnetic mass 
difference for pions and kaons in [19] and related work can be found in [6, 20, 21]. We 
obtained very nice matching for four-point functions in the presence of quark masses. 
At present, this non-trivial matching has only been obtained in the 1/A^c-technique. 

A point we do not discuss here is the choice of gauge for the X-boson. For 
reasonable gauges the effect are suppressed by extra powers of 1/ M"^ so they arc not 
needed here. In any case the discussion for the photon in [19] can be easily extended 
to show that there is no problem here in general either. 

We want to calculate some matrix element between hadronic states in general. 
For our example in (3.4): 



^ _9a^(K0\ I I d'y^L^{x)L,{y)Px{x,y)\K' > . (3.22) 

The equality follows from the discussion in the previous subsection. Px{x,y) is the 
X-boson propagator 

^'x(...)^/gi^;^e-<-'«. (3.23) 

The matrix element corresponds to a four-quark operator in an OPE in the X-boson 
effective theory. 

We now rotate the integral in (3.23) into the Euclidean and split the integral 
into two parts \pe\ > IJ' and \pe\ < IJ^- One of the reasons to rotate to Euclidean 
space is that in the lower part then all components of are also small. Also, in 
general, amplitudes are much smoother in the Euclidean region, so approximations 
are generally better behaved since threshold effects and the hke are smeared out. In 
particular we split 

I d'^pE = I i [d\pE\ \pe? + Pd\pE\ \peA . (3.24) 
long-distance short-distance 
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This is the same procedure as apphed to the electromagnetic mass difference of pions 
and kaons in [19]. We choose Mx ^ ^ such that in the long-distance part we can 
neglect all momentum dependence of the X-propagator and we expand the short- 
distance part in 1/M^. 

Form-factors of mesons at large Euclidean momenta p\ are suppressed by 
so in the short-distance part the high momentum has to flow back through quarks 
and gluons to leading order in 1///^. The short-distance part is a two-step calculation, 
we evaluate the the integral with the momentum flowing back through gluons and 
quarks using the box diagrams of Figure 3 only. The vertex diagrams do not involve 
a large X-boson momentum so they are part of the long-distance calculation. Here 
there is no infrared ambiguity, the possible infrared divergence is regulated by /x, the 
lower limit of the short-distance integral. The result is: 

mK')sD-x = ^^^mS. + 52|i^°)^7i In ^ . (3.25) 

Here as is already suppressed by 1/Nc so it is sufficient to calculate {K^\Si + S2\K'^) 
to leading order in l/Nc, i.e. oo. 

Now the long-distance part can be calculated also order by order in l/Nc- The 
leading order in 1/Nc is 

(^l^°)LZ.-x-leading = + ^2|i^°>iv.-.oo (3.26) 

where the subscript Nc ~* oo indicates again the leading in 1/Nc contribution. For 
the subleading in 1/Nc part it is sufficient to replace l/{p^ — M^) by —l/M^ so we 
obtain 



2Mx 



{Tr\Jd^xJd'yAL^{x)L,{y)PE{x,y)\K''>^,^^^^^ a.ii,) ^ /i 



(3.27) 

The subscript 1/Nc sup means only the 1/Nc suppressed part and 

PE{x,y) = £d\pE\ \PE? j df)p,e-*(^-^)-^'-. (3.28) 

We can now insert the value of g\s=2 of Eq. (3.21) into (3.27) and we see that 
all Mx dependence disappears. The final result is 

+ m j dS j d'y4L^{x)L^{y)PE{x,y)\K' >>tf^^^^^ \ (3.29) 
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The V and ^lq dependence had aheady disappeared in the previous section. The 
/X dependence that remains in the overall coefficient and in the upper limit of the in- 
tegral in Pe{x, y) are both at the same order in l/N^ and will cancel if the low-energy 
approximation to QCD used to evaluate the long-distance l/Nc suppressed contri- 
bution matches onto QCD sufficiently well. The complete short-distance corrections 
to two- loops can be seen in Figure 2. 

Notice that this proof only required the use of the 1/A^c-expansion at a rather 
late stage, thus the arguments all go through provided the scales /i, /ic, and u are 
chosen such that no large logarithms of their ratios appear. The identification of the 
short-distance pcrturbativc QCD quantities with the quantities in the long-distance 
1/Nc suppressed part need only be done at the level of two-quark currents. Once 
this is done, and as mentioned before this is normally a requirement for them, the 
scheme dependence introduced in short-distances is fully eliminated. 

The only, but difficult, obstacle for an almost perfect matching is now to find a 
model that is viable up to scales fi where perturbative QCD is applicable. 

The matching between long- and short-distance can be systematically improved 
with low-energy realizations of QCD which are better and better at high energies. 



3.3 Gluon Gauge Dependence 

In the above we have always used the Feynman gauge. So what happens now in other 
gauges. We can investigate the gauge-dependence of both matchings done before. 

The long-distance-short-distance matching in Section 3.2 can be easily shown to 
be gauge-independent to the order we are working. 

The transition to the X-boson scheme from the operator product expansion is in 
fact dependent on the gauge. Actual calculations show that the gauge dependence is 
purely a long-distance phenomenon^. The problem is that the state we have chosen 
in Eq. (3.6) is not a physical state. We could have chosen an infrared well-defined 
observable to do the matching from the OPE to the X-boson theory. This would 
mean going to on-shell massless quarks and including soft-gluon radiation as well. 
Then the gauge independence would have been explicit. 

The other alternative is to use a gauge-independent infrared regulator but with 
it fixed at the same value for both sides or simply, as we did, to use the same gauge 
and infrared regulator on both sides of Eq. (3.6). We have checked that the gauge 
dependence with our off-shell quarks is the same on both sides and it does cancel in 
the value of g\s=2- 



^See also the short discussion in the second reference in [9], Section III.C. 
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4. Scheme Independent Results for the Bk Parameter 



The kaon bag-parameter Bk is defined by 



1 

{K^\K') ^ -iC^s=2^fWKBK . (4.1) 



The fully scheme-independent result from (3.29) gives 



7i In ( - ) + ri - fi 



BM ; (4.2) 



where 



with 



„ , , 3/ , {K«\li'xld'yiL-(x)L-{y)PE(x,y)\K''>t,^^„; 



In our previous work [12] we used, as mentioned there, fi = and r^^^ — —7/6 [10]. 
Here we adjust for that. 

The remaining matrix element can now be calculated in various ways. If we use 
NLO Chiral Perturbation Theory we obtain in the chiral limit^ 

Away from the chiral limit the expression are much more cumbersome but for the 
lowest order CHPT they can be found in [12]. This way we can obtain it as a series 
in n'^. But as plotted in Figure 4 we see that this series is already breaking down at 
values of /i ~ 500 MeV. We therefore need an extension of CHPT that includes the 
above result as a low energy limit. The ENJL model has the same chiral structure as 
QCD and also includes spontaneous chiral symmetry breaking. It thus has CHPT as 
an automatic built in limit. We then choose the three ENJL parameters such that the 
CHPT parameters to order are well described so Eq. (4.5) is included. This version 
of the model also correctly reproduces a lot more hadronic phenomenology [23] . In 
Section 5 we summarize some of its advantages and disadvantages and how we expect 
the latter to be of little influence for the numerical result. The results for Bj^{n), the 
chiral limit result, and Bk{^)^ the result at the quark masses corresponding to the 
physical pion and kaon masses arc given in Table 1 of [12]. We have plotted these as 
well in Figure 4. We can see here very well that the ENJL model includes the CHPT 
results at low and considerably improves on it at high /x, both in the chiral limit 
and for physical values of the quark masses. All curves are in the nonet limit. 

^The diagrams contributing are those of Figure 3 in [12] but the vertices can now also be those 
from the order CHPT Lagrangian. 
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Figure 4: The long-distance part of the X-boson exchange calculation. We show Bk{ij) 
of Eq. (4.2). The lines are leading- (LO) and next-to-leading-order (NLO) CHPT and the 
ENJL model for the chiral limit and LO CHPT and the ENJL model for the physical quark 
masses. Notice the improvement of the ENJL model over CHPT. 

We can now put these together with the short distance part with ^^(M^) = 0.334 
and 11 = v. We use for definiteness the NDR results for gAs=2 but we could have used 
the HV ones with the same answer, see Figure 2. We also included for the massive 
case the needed correction of 0.09±0.03 for the octet versus nonet case[8, 12], mainly 
due to rj-r]' mixing. The result is shown as a function of u in Figure 5. Notice the 
quality of the matching. 

In [17] a better identification of the scale together with the Leading-Order CHPT 
approximation was used, this corresponds to keeping the quadratic divergence only. 
In the chiral limit this correspond to the first two terms in (4.5). Negative values 
for Bk and the related 27-plet coupling when using this approximation in the chiral 
limit (see Figure 5), were also seen there[17]. Notice in the same figure how much 
the use of higher orders estimated via the ENJL model improves the matching in 
the chiral limit and outside the chiral limit [12]. The variation is more due to the /i 
dependence since as shown in Figure 2 there is little dependence on u. The massive 
case has an extremely good stability for 700 MeV to 1 GeV while the chiral case is 
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stable for 550 to 700 MeV. 



LOCHPTchiral 
ENJLchiral 

LO CHPT mass 
ENJL mass 




Figure 5: The results for Bk in the NDR scheme with asiMr) = 0.334 as a function of 
fj, = v. The curves are from top to bottom: ENJL massive, leading-order CHPT massive, 
ENJL chiral limit, leading order CHPT chiral limit. Notice the ENJL stability as compared 
with leading-order CHPT. Next-to-leading-order CHPT is even more unstable. 

As already noted several times before, [8, 12] and references therein, the non-zero 
quark-mass corrections to Bk are quite sizable. 

We can now study the dependence on the variation with the input. In Figure 
6 we have plotted the result with the same input for the HV scheme with the same 
value for asiM^-) and in the NDR scheme for ^^(Mt-) = 0.36 and 0.31, including the 
value of [24]. We took again ^ = v. The stability is essentially unchanged while the 
actual values are very similar. 

Finally we want to check the variation with [i ^ v. This again provides a check 
on neglected corrections of order o?^. In Figure 7 we plotted the NDR case for 
Q;s(M.r) = 0.334 as a function of \i for, from top to bottom, v = (1.5, 1.2, l,0.8)/i. 
The last curve we did not plot for low values of fi since this gets very large values 
for Again, the variation is rather small. 

The uncertainty can be judged from all these variations. We use the variation 
with asiM-r) and add a similar error for the remaining model uncertainty to obtain 



14 



1 I 1 1 1 1 1 1 r 



< m 




Q I Li I I I I I I I I I 

0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

V (GeV) 



Figure 6: The variation of Bx with as as a function of u. The top set of curves is the 
massive case, the bottom set is the chiral hmit case. Curves are from top to bottom: NDR 
with asiMr) = 0.31; HV with a,(M^) = 0.334; NDR with a,(M.^) = 0.334; NDR with 
asiMr) = 0.36. 

the result given in Eq. (6.1). 

5. Some Remarks on the ENJL Model and Possible Improve- 
ments 

The model we use has the chiral structure of QCD at large N^, i.e. it is a left-right 
symmetric model which breaks down to the vector subgroup. Quark masses break 
this symmetry as in QCD. 

It provides a picture of build-up of constituent quarks out of massless ones. It 
has two drawbacks, it doesn't confine and in some channels it does not reproduce the 
same high-energy behaviour as QCD. These violations are small at low-energies. It 
reproduces CHPT to order [23] , and provides an estimate of all the higher order 
corrections. It reproduces a large amount of low-energy phenomenology for both the 
anomalous and the non-anomalous sectors [23] with an accuracy at the level of 20% 
to 30%. 
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\i (GeV) 

Figure 7: The variation of Bx in the NDR scheme for, from top to bottom, u = 
(1.5, 1.2, 1, 0.8)//. The top set of curves shows the massive case, the bottom set the chiral 
hmit. 

The drawbacks we think will not affect our numerical results too much because: 

1. We do all our integrations in the Euclidean domain. Here the effect of sharp 
states are smeared out so a model without confinement that reproduces smeared 
quantities correctly should be all right. 

2. We only use the ENJL model at intermediate momenta where CHPT fails and 
QCD might not work yet. The good matching we obtain is an indication that 
in the regime where we use it the ENJL model works satisfactorily. 

The model works all right but is still the weak link in our approach. Work on 
some matrix-elements using more direct QCD based arguments has been done [21, 
25] and alternatively, we can augment the ENJL model in ways that improve the 
matching with QCD by requiring matching in as many relevant Green functions as 
possible. 
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6. Conclusions 

In this paper we have shown how the short-distance large logarithms can be re- 
summed in a consistent fashion while retaining the beauty of the approach. We 
have thus shown how to realize Bardeen's picture of tracking scheme-dependence all 
the way [4, 5] from the W-boson mass to very low scales. The approach is general 
but we illustrated it in the case of the Bx parameter explicitly. 

The main assumption is that we know how to hadronize quark currents, once 
that is done, our approach tells how the four-quark operators should be hadronized. 
This can be systematically improved by higher order calculations within the pertur- 
bative domain. The low-energy hadronic realization of currents of QCD is the only 
remaining model dependence. At very low energies we know that CHPT describes 
the QCD behaviour correctly, thus this should be included. In the intermediate do- 
main we can try to use data as much as possible and/or models with various QCD 
constraints. We have argued that the ENJL model is a good first step beyond CHPT 
and our final results show this improvement. 

By varying a^, the matching scales and conditions, i.e. /j, ^ u we can get an 
estimate for the error involved. Our final result for is: 

Bk = 0.77 ± 0.05 (as) ± 0.05 (model); 

B^ = 0.32 ± 0.06 (as) ± 0.12 (model). (6.1) 

The difference between these results and the ones in [8, 12] is only the short-distance 
scheme dependence. Our results here are short-distance scheme independent. The 
main uncertainty are the value of ^^(Mr), and a similar error for the remaining 
model dependence. The small model dependence error we quote is due to the almost 
cancellation of the l/N^. corrections, remember B^ = 0.75 at leading order in l/N^.. 
Possible systematic errors are of course difficult to estimate. The scale and scheme 
dependence are consistently matched at order for current x current operators and 
are only a small part of the final error, especially for the physically relevant massive 
case. 

The scheme dependence has not been discussed in other non-leptonic weak matrix 
elements calculations like the 1/Nc technique in [17] or the chiral quark model in [26]. 

We will apply the same procedure to the AS = 1 transitions where first results 
on the A7 — 1/2 rule were obtained in [12] using r^^^ from [10] and fi = 0. Work 
is in progress to perform the remaining calculations also for AS" = 1 transitions and 
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